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We present a control-theoretic analysis of the system consisting of a two-level atom coupled with 
a quantum harmonic oscillator. We show that by applying external fields with just two resonant 
frequencies, any desired unitary operator can be generated. 
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I. INTRODUCTION 

In this paper, we apply theoretical concepts of quan- 
tum control to the joint system consisting of a two 
state system coupled with a quantum harmonic oscilla- 
tor. Such systems are ubiquitous in Nature. For example, 
coupled atom-oscillator systems form the basis for the 
ion trap quantum computer [l[. Other examples include 
a single atom in a cavity Q, a super-conducting qubit in 
a cavity 3, and control of single atom lasers (HQ. In 
Law and Eberly showed that arbitrary states can be syn- 
thesized by using just two resonant frequencies, a result 
experimentally verified in (23|, and [TtJ showed that the 
two-level atom-oscillator system could be controlled by 
fine-tuning the Lamb-Dicke parameter. Here we prove 
that the dynamics of such systems is controllable with- 
out any fine-tuning or special state preparation: with 
the proper sequence of pulses, it is possible to perform 
any desired unitary transformation on the Hilbert space 
spanned by the atomic states together with the lowest n 
energy levels of the oscillator. 

In this paper, we will use the ion trap as our model 
system. An ion trap quantum computer can be mod- 
eled as a collection of N particles with spin | in a one- 
dimensional harmonic potential. Laser pulses incident 
on the ions can be tuned to simultaneously cause in- 
ternal spin transitions and vibrational (phonon) excita- 
tions, thus allowing local internal states to be mapped 
into shared phonon states. The computational qubits are 
encoded by two internal states of each ion and the collec- 
tive vibration of the trapped ions acts as the information 
bus. In this manner, quantum information can be com- 
municated between any pair of ions and logic gates can 
be performed. Several key features of the original pro- 
posal in [l[ , including the production of entangled states 
and the implementation of quantum controlled opera- 
tions between a pair of trapped ions, have already been 
experimentally demonstrated (see, e.g., d, d, El EH)- 
Meanwhile, several alternative theoretical schemes (see, 
e.g., [H El EI EE El E3) have also been developed 
for overcoming various difficulties in realizing a practi- 
cal ion-trap quantum information processor. All these 
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proposals either require fine-tuning of the Lamb-Dicke 
parameter or an initial cigcnstatc of the vibration mo- 
tion. Here we present a control theoretical analysis and 
show that in the Lamb-Dicke regime by using two res- 
onant frequencies, any unitary transformation within a 
finite level of the harmonic oscillator can be generated. 
Unlike, e.g., [13], no fine-tuning of the Lamb-Dicke pa- 
rameter is required to obtain complete control. While the 
proof of controllability is somewhat involved, because of 
the fundamental nature of the system to be controlled 
and because of the wide range of potential application, 
we present this proof in detail. As will be seen below, 
the difficulty of the proof arises because, in the absence 
of controllability of the Lamb-Dicke parameter, one must 
combine discrete and continuous control theoretic tech- 
niques. The current proof can be regarded as extending 
the techniques of the paper 0, E3| from controlling 4 
states to controlling m states, where m can be arbitrary 
large. We begin in section |TT] by presenting the usual 
Jaynes-Cumming model for spin boson interaction. Wc 
then make the controllability analysis of the system in 
section Mil 

II. LASER-ION INTERACTION MODEL 

The physical situation we consider is a two-state atom 
(frequency lu c ) coupled to a harmonic oscillator (tran- 
sition frequency u z ), driven additionally by an exter- 
nal field (frequency u>). We will follow the ion trap 
model @, EH]- The free Hamiltonian of this system is 
Hq = fiui c ^j- + hu> z a^ a, where a z is a Pauli spin operator 
and a annihilates a phonon. Turning on the electromag- 
netic field of a laser gives an interaction Hamiltonian 

h! = -ji-b, (i) 

where fl = [mj/2 is the magnetic moment of the ion and 
B = Bx cos(kz— cut+Q) is the magnetic field produced by 
the laser. Here z = zo(a + a'), where Zq = \Jh/2Nmui z is 
a characteristic length scale for the motional wave func- 
tions and m is the mass of an ion. 

We consider the regime in which rj = kzo <C I . In this 
regime, we may determine the effect of a laser pulse at a 
specific frequency u) by expanding Eq. {T]) in powers of rj 
and neglecting rapidly rotating terms. Then pulsing on 
resonance (ui = oj c ) allows one to perform the transfer- 
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mat ion 



R{9, (/)) = exp [t0(e**<j + + e _i *0] 



and pulsing at the red sideband frequency (u = io c 
gives 
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exp [^(e^cr+a 1 
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(2) 
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In each case, the parameter 6 depends on the strength 
and duration of the pulse and <fi depends on its phase. 
In the next section, we show that by just using these 
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FIG. 1: Graphical representation of the quantum harmonic 
oscillator driven by a sinusoidal resonant field fields uo c and 
uj r — uic — as shown. The strengths of the u) c transition 
couplings are independent of the harmonic oscillator quantum 
number n, whereas the strength of the u) r transition couplings 
increase as the square root of the quantum number n. 

two frequencies any unitary operator can be generated. 
The basic idea in proving controllability is an extension of 
[lH , [lj| . Using the feature that the transition frequencies 
increase as the square root of the quantum number, we 
apply only pulses that leave the system confined within 
the Hilbert space spanned by the first n oscillator lev- 
els. This requirement means that the only a discrete set 
of pulses can be applied at the red sideband frequency. 
Meanwhile, a continuous set of pulses can be applied at 
the resonance frequency. As a result of the use of both 
discrete and continuous controls, the resulting control 
problem is technically somewhat involved. Nonetheless, 
it can be solved completely, as we now show. 



III. CONTROLLABILITY ANALYSIS 



We denote E pq be the matrix such that has all the 
entries equal to zero except the pq entry, which equals 1. 



It is easy to check that E pq E rs 



— SgE pS . 



The Hamiltonian, after absorbing the imaginary num- 
ber i, can be represented as skew-Hcrmitian matrices. If 
we take the eigen-state of the free Hamiltonian as the ba- 
sis, then after re-scaling the time unit, the various Hamil- 
tonians in the interaction frame can be represented as 
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When the Hamiltonian H3 or H4 is applied, | t)l TO ) is 
connected to | l)\m + 1). We restrict the evolution time 
T under these two Hamiltonian to satisfy T^/m = kir, 
while k is integer, so the subspace of states spanned by 
{| l)T)b')|j < m } is preserved. We show that under 
these restrictions any unitary matrix within any finite 
harmonic level still can be generated. 
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SU(4) 



Choose k = 2p + 1, we can have 



Let's first work out the case of m = 1, show that we 
can generate SU (4) on the subspace spanned by states 

|l)|0},|T)|0),||)|i},|T)|i) 

Restricted to this subspace, 
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and the unitary operators we can generate using H3 and 
H4 are 
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Ui = R 
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Since 17(4) is compact, the infinite sequence 
{t/i, XJ\ , £7f , C^i--.} has a convergent subsequence, 
i.e., there exists p\ > p 2 <G N such that J/f 1 — £/f 2 is 
arbitrary close to zero, when this is true then jjP 1 ^ P2 ^ 1 
is arbitrary close to ?7f , i.e., ?7f can be approximately 
generated to arbitrary accuracy. But 
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choose k such that is arbitrary close to an integer, we 
can get the transformation 
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Choose k — 2p, by varying p, R (^7=, § ) forms a dense 
subset of the one parameter group 
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Thus we have the generator 
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, add it to H 2 , 
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Subtracting this from Hi and dividing by a factor 2 yields 



Ha 
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Similarly by using H 2 and U\, we can get 



H 7 
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We now show that {H 5 , H 6 , H7} generate all the skew- 
Hcrmitian matrices on the subspace. First 



N-2 



H$ — H 5 - H 7 — E, 



k(k+i) - E( k+ i) k 



k=i 



Here N = 4, we will do the following computation using 
the general N, as this will be used for the proof of the 
general case. Now, 
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We first show that H-j and Hg, generate all the real skew- 
symmetric matrices of size N x N [2(| , let 

M N _i = H 7 = S(jV-l)AT — ^AT(W-l) 
Mn-2 =Ws,M N -i] = E( N _ 2 )N — E N ( N _ 2 ) 

=[H&i + — E( N _ 3 -j N — -Siv(jv-3) 

A/jv-4 =[-#8, Mat_ 3 ] + Mjv-2 = S(JV-4)W ~ ^JV(/V-4) 

Mi =[iT 8 , M a ] + M 3 = £ijv - £wi 

and [Af„M g ] = E qp - E pq , Vp ^ q e {1,2,..., N- 1}. 
Thus we can generate complete basis for skew-symmetric 
matrices. Similarly 

Jn-i = H 6 = i(-E(Af_i)Ar + E N ( N _i-j) 

Jn-2 =[Hg, JjV-l] = i{E(N-2)N + -EW(iV-2)) 

Jns =[Hg, Jn-2] + Jn-i = i(E(N-3)N + £jv(iV-3)) 



Ji =[H&, J2] + J3 — i{E\N + Ejsri) 

and 

[M q ,J p ] =i(E qp + E pq ) 

[^(,E qp -\- E pq }, E qp Epq\ 2%{E pp E qq *) 

\/p 7^ q G {1, 2, ...,N— 1}. So we can generate a full basis 
for all N x N skew-Hcrmitian matrices. This proves the 
S77(4) case. 

B. General case 

Now we generalize our proof to the controllability on 
SU(n) for any n. It is not necessary to check the case for 
each n, as SU(n{) is a subgroup of SU(ri2), for m < n,2, 
the controllability on SU(n 2 ) implies controllability on 
SU(n\). It is sufficient to prove the result for infinitely 
many rij as n, — > 00. 

Take the subspace up to Harmonic level to, i.e., 

{U)|0),|T>|0>,U>|l),|T)|l),...,U)|m),|T)|m)} 

where (m — 1, m+ 1) are both prime. We shall prove the 
controllability on this subspace. The twin prime con- 
jecture claims there exists infinitely many such primes. 
If the twin prime conjecture is false, then the following 
proof works only up to n = 2m + 2, where to is the 
largest known twin prime. As of 2006, the largest known 
twin prime is 1003 145 1 25 4 40 1 5 -2 171960 ±1, which is large 
enough for most physical systems. Below, we generalize 
the twin-prime proof to show controllability for all n. 

If we restrict the evolution time T for H4 to satisfy 
T\Jm + 1 = kir, where k is an integer, then the angle ro- 
tated between | T)b — 1) an d | l)\p) is ^JpT = kJ-^jir. 



We divide the numbers {1,2, ...,m} into groups Gi,i = 
1,2,..., such that in same group Gi the angles rotated 
under the above evolution are rationally related to each 

other, i.e., pi, P2 are in same group, if and only if is 

a rational number. For example, {1, 1 • 2 2 , 1 • 3 2 , 1 -p 2 } 
forms a group, where p\ < to, (pi + l) 2 > to, similarly 
other groups are {2, 2 • 2 2 , 2 • 3 2 , 2 • pf},{3, 3 • 2 2 , 3 • 
3 2 , 3 • p 2 }..., specially {to — 1} itself forms a group. 

As to + 1 is a prime number, k y / ' ^^(mod2) are irra- 
tional numbers for all p < to. Accordingly, we can vary 
k such that, except the angles relate to one group Gi, all 
the other angles are arbitrary close to zero. This way we 
can construct the generator 

H * = ^2 Vj( E 2j(2j+l) - E {2j+1 )2j) 

Add all Hi to H 2 , we get a matrix similar to H 5 in the 
SU(4) section, with only nonzero entries at the first off- 
diagonal. Denote this matrix by H$. 

To prove the controllability, we just need to show that 
we can also generate matrices similar to H§ and H7, i.e., 

E(N-1)N — En(N-1) 

and 

i(E(N-l)N + ^AT(JV-l)) 

here N — 2m + 2. 

As {m— 1} itself forms a group, say Gj , we can generate 

Si = / 1 Hj = E(2 m — 2)(2m-l) ~ ^(2m-l)(2m-2) 

bracket it with H 2 = Ylk=o E (2k+i)(2k+2) -E {2k+2){2k+1) , 
we can get 

52 =[H 2 , S\] 

=-B(2ra-3)(2m-l) — -E-(2m-l)(2m-3) 
+ -E'2m(2m-2) — -E(2m-2)2m 

then bracket S 2 with S\, 

53 = [Si, S 2 ] 

— S(2m-3)(2m-2) — "(2m-2)(2m-3) 
+ S(2 m -l)2m — ^2m(2m-l) 

we sec that S3 is nothing but the restriction of H 2 on the 
subspace spanned by 

U)|to-2),|T)|to-2),|4)|to-1),|T}|to-1) 

Similarly Hi can also be restricted to this subspace. 
From the SU (4) case, we know we can generate any skew- 
Hermitian matrix on this subspace, specifically we can 

have S4 = £7(2m-l)2m ~ ^2m(2m-l) 
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Now pick the group G p to which m belongs. We get 

Hp = X! Vj( E 2](2j+l) ~ E( 2 j+l)2j) 

=-E'2m(2m+l) — ^(2m+l)2m (12) 
+ V7(-^2j(2j+l) - #(2j+l)2j) 

Bracket S4 with since all the numbers in G p have the 
form the second term in the right side of the above 
equation commute with S4. Accordingly we obtain 

S 5 = [Si, H p ] = B(2 m -l)(2m+l) — ^(2m+l)(2m-l) 

Now, bracket S5 with S4, 

Se — [85,84] — E2 m {2m+1) — ^(2m+l)2m 

comparing with S\ , we see that we just moved one block 
down. Repeat what we did with Si to Sq, we can get 

87 = -E-(2m+l)(2m+2) — ^(2m+2)(2m-l) 
= ^(JV-l)JV — £jV(JV-l) 

This is the matrix we need to generalize our proof of 
controllability on SU(4). Similarly, we can get 

i(E(N-l)N + ^JV(JV-l))) 

together with i/ 5 , we are able to generate all the skew- 
Hcrmitian matrices of size N x N, which proves the con- 
trollability on SU(N). This completes the proof: driving 
the fundamental frequency and the red sideband suffice 
to control the two-level atom coupled to an harmonic os- 
cillator. 



Remark 1 From the proof we see that the only two 
properties of the pair (m — l,m + 1) we used are: 

1: \J^fi arc irrational for all p < m. 

2: There exists one group consists of only one number. 
It is convenient to pick twin primes, but there exist other 
choices. For example, we can choose m+ 1 = 2q, where q 
is an odd prime. Under this choice condition 1 still holds, 
and q itself forms a group. So our proof, while expressed 
in terms of the twin prime conjecture, actually holds for 
all n. 



IV. DISCUSSION 

We have proved the controllability of the dynamics of 
the coupled two-level system/harmonic oscillator. Be- 
cause of the discrete nature of the controls, the proof 
was somewhat involved. In addition, the system is only 
fully controllable in the limit that the number of control 
pulses goes to infinity. In any realistic setting we will 
have only a finite time and a finite number of pulses that 
we can apply. The question of the rate of convergence 
of such discrete schemes is an important open question 
in control theory and in quantum information, and will 
be investigated elsewhere. For the moment, we note only 
that accurately generating arbitary members of SU(4) 
and SU(n) for n < 10 or so via the techniques described 
here is well within the reach of current experiment. 
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